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Abst  ract 

Integral  expressions  for  the  pump  and  generated  fields  are  presented 
here  for  the  case  of  second  harmonic  generation  of  a  focused  Gaussian 
pump  beam  incident  on  a  nonlinear  crystal  .  The  birefringent  walk-off  of 
the  generated  beam  and  the  effect  of  pump  depletion  are  included  in  the 
theory. 


1  Introduction 

The  theory  of  second  harmonic  generation  (SHG)  of  a  focused  Gaussian  beam 
in  the  presence  of  beam  walk-off  presented  by  Boyd  and  Kleinmaid1)  has  been 
widely  used  in  the  nonlinear  optics  community.  A  limitation  of  their  method  is 
the  exclusion  of  pump  depletion  which  makes  it  inapplicable  to  the  case  of  high 
conversion  efficiency.  The  numerical  method  developed  by  Smith  ^  and  recent 
work  by  Wang  and  Weiner  fb  and  by  Kasamatsu,  Kubomura  and  Kan  ■ 1 ■'  ex¬ 
tend  the  theory  of  Boyd  and  KJeinman  to  include  the  effect  of  pump  depletion* 
We  present  here  an  alternate  technique  to  calculate  the  spatial  distribution  of 
the  pump  and  the  SHG  fields  in  presence  of  diffraction,  linear  absorption,  phase 
mismatch,  beam  walk-off  and  pump  depletion.  The  field  distributions  are  ex¬ 
pressed  as  multiple  integrals,  which  are  considerably  simplified  for  the  case  of 
collimated  beams,  i.e.,  when  the  spread  of  the  pump  beam  due  to  diffraction 
inside  the  crystal  can  be  ignored.  This  is  often  the  case  when  high  peak  power 
pulsed  lasers  are  used  as  pump  beams  because  the  pump  is  intentionally  kept 
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collimated  through  the  nonlinear  optical  medium  to  avoid  damage  to  the  mate¬ 
rial. 

In  the  theory  presented  here  the  incident  pump  beam  is  assumed  to  be  fo¬ 
cused  at  the  crystal  center  -  this  restriction  can  be  easily  removed  later.  The 
aim  of  this  work  is  to  describe  the  general  approach  taken  and  to  provide  ex- 
prcssions  for  the  electric  holds  expanded  up  to  fifth  power  in  the  second-order 
nonlinear  optical  coefficient,  d.  The  method  described  here  is  an  extension  of 
an  early  work  on  optical  parametric  oscillation  (OFO)^5)  and  can  be  readily 
applied  to  other  nonlinear  optical  processes  of  current  interest,  such  as  pump 
resonant  SRO  or  resonant  second  harmonic  generation  (SHG).  As  pointed  out 
in  ref. 3,  the  Green’s  function  method  is  easily  extended  to  the  time  dependent 
case  as  well.  However,  a  limitation  of  the  perturbative  approach  used  here  is 
that  multiple  integrals  need  to  be  evaluated,  which  becomes  computationally 
time  consuming  when  six,  eight  or  higher  dimensional  integrals  are  involved. 

In  the  following  sections,  the  derivation  of  the  relevant  wave  equations  is 
provided  first,  to  set  up  the  starting  equations.  Then,  the  Green’s  function 
method  of  solution  in  terms  of  a  power  series  in  the  nonlinear  coefficient  is 
described,  and  the  expressions  for  the  pump  and  the  SMG  beam  electric  fields 
are  evaluated  for  a  fewT  terms. 


2  The  Starting  Equations: 

The  notations  for  electric  fields  and  nonlinear  polarizations  used  here  will  gen¬ 
erally  follow  the  convention  adopted  in  the  Handbook  of  Nonlinear  Optics ^ 
(using  the  SI  system  of  units).  The  subscripts  p  and  s  denote  the  pump  and  the 
SIIG  fields,  respectively.  The  fields  are  assumed  to  be  linearly  polarized,  and 
only  the  scalar  components  of  the  electric  fields  are  considered  in  the  nonlin¬ 
ear  interaction.  The  equation  16.4.5  in  Quantum  Electronics ^  is  taken  as  the 
starting  point  of  the  analysis  here: 

V  E  =  ,^— +w— +  ,i0— PNL  (1) 

The  electric  fields  at  the  pump  and  the  SHG  frequencies  {up  and  u/Si  respec¬ 
tively)  are  written  as 

£J“,’(r1 1)  =  Ep(r)e~iu”t  +  cc  (2) 

£“'■  (r,  t)  =  Es(  +  cc  (3) 

arid  the  total  nonlinear  polarization  as 

PnUt,  t)  =  (P&tte-**  +  cc)  +  +  cc)  (4) 

Following  ref.  6  (chapter  2)  we  assume  that  the  nonlinear  polarization  terms  at 
the  pump  and  the  SHG  frequencies  can  be  written  as 

V?fL(r)  =  4di0Ea(r)E;(r)  (5) 
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T%*l(t)  =  2<k0Ep(r)2  (6) 

Inserting  equations  2.3 ,4,5  and  6  in  equation  I,  and  collecting  terms  oscillating 
at  frequencies  ujp  and  ,  we  obtain. 

(V2  +  k2pi)Ep (r)  -  -*,w*7#L(r)  (7) 

(V2+^)£!(r}  =  -/io^;iW  W 

where 

kpi  —  k2  +  Wp/ifjcrp  -  k2  +  iwaju0crs  (9) 

In  the  above,  it  has  been  assumed  that  in  equation  1,  tp8  —  ton2,s  and  that  the 
wavevectors  are  kp  s  —  71p^ojp,3/c  where  np  and  n8  denote  the  refractive  indices 
of  the  medium  at  the  two  frequencies. 


3  Green’s  Function  Solution 


The  solution  to  the  inhomogeneous  differential  equations  7  and  8  above  is  ob¬ 
tained  using  the  Green's  function  method: 


Ep(r)  =  - 


Es(r)  =  - 


^0^ 

4tt 


4tt 


£/rfr'V^(r) 

R„ 


(10) 

(11) 


where  Rp  and  R8  denote  the  magnitudes  of  the  displacement  vectors  connecting 
the  points  F  and  r  for  the  pump  and  the  SHG  fields,  respectively.  The  pump 
beam  is  assumed  to  propagate  in  the  medium  as  an  ordinary  wave,  so  that 


RP  =  I  r  -  r' 

=  y/(x  -  z')2  +  (y  -  +  (z-  (12) 

If  the  SIIG  field  propagates  through  the  nonlinear  crystal  as  an  extraordinary 
wave,  say  with  a  walk-off  angle  denoted  by  /),  and  if  the  Cartesian  coordinates 
are  chosen  such  that  the  propagation  vector  of  the  SHG  field  lies  in  the  x,  z 
plane  and  is  at  angle  p  with  the  z  axis,  we  have  ^ 

R»  =  V{(^  -  P*)  -  {%'  -  Pzf)Y  +  {V  “  V1)2  +  iz  ~  z')2  (13) 

For  small  absorption,  equation  11  can  be  re-written  as: 

kpi  ^  kp+i ^  k31  S  ks  +  ty  (14) 
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where  aPrft  =  crypto /ePtS  denote  the  the  linear  absorption  coefficients*  The 
Fresnel  approximation 

Rp-(z-  z')  +  j(rrjrj{(x  -  x')2  +  {y-  v')2}  (is) 

and 

R*  =  (*  -  *0  +  2(7T7)  «*■  -  *'?  +  is  -  v')2}  («) 

is  used  in  the  exponents,  with  xa  —  x-p(z-zi)  and  the  approximation  Rp ,  Hs  = 
{z  —  z*)  is  used  in  the  denominators.  Equations  10  and  11  then  reduce  to 


,  _  fto^p  f  dx!dyfdz{ 

4  7T  J  Z  -  Zf 

eikp(z-z/)-^{z-zf)+~£^{(x~x,)2+(y-yf)2}<pUp 


(17) 


and 


EM 


f 

4n  J 


dxrdyfdz/ 

z-zf 


ik,(x-Zr)-2f(X-z')+ 


{(i„-a:')2+(v-^)2} 


(18) 


The  incident  pump  beam  is  assumed  to  have  Gaussian  cross-sections  and  to  be 
focused  at  the  center  of  a  nonlinear  medium  of  length  l.  Assuming  the  origin  of 
the  coordinate  system  to  be  at  the  incident  face  of  the  medium,  the  amplitude 
Epli}  (r)  of  the  incident  pump  field  can  he  written  as 

—  _ ^P°  cik„z-g„z/2c-(x2+y2)/{wZr(l+2i(z-if)/b„)  ^qgj 

p  1  +  2 i(z  -  | )/bp 

where  bp  =  kpw^p  denotes  the  confocal  parameter.  In  the  first  order  approxi¬ 
mation,  he.,  when  pump  depletion  is  ignored,  it  is  assumed  that  the  pump  field 
Ep  remains  equal  to  Ep°\  a  nonlinear  polarization  at  frequency  w*  is  generated 
through  equation  6,  which  gives  rises  to  an  electric  field  (dEj1^)  at  frequency  usB 
through  equation  1L  This  generated  SHG  field  interacts  with  the  pump  field 
Ep  to  generate  a  nonlinear  polarization  through  equation  5,  which  generates 
(through  equation  10)  a  component  of  the  electric  field  at  the  pump  frequency 
(denoted  by,  say,  d2E^).  Following  this  procedure,  a  series  of  SHG  and  pump 
field  terms,  respectively  odd  and  even  powered  in  d,  can  be  obtained.  That  is, 
in  general,  the  SHG  and  the  pump  fields  can  be  written  as 


Es  -  dE '»  4-  d3E<3)  +  dbEib) . . . 

(20) 

Ep  -  Ep  +  d2Ep  +  diEp  +  dti  Ep  . . . 

(21) 

Substituting  these  series  in  the  equations  5  and  6,  the  nonlinear  polarization 
terms  can  also  be  expanded  in  power  series  in  d : 


where 


and 


sr-JLlJf  sTJfi  ^  U  1  sTJ  I  ^  ^  ^ 

yNL  -  hnl  +  rNL  +  rNL  .  - . 

(22) 

“ NL  —  ^ NL  ^  rNL  +  rNL  ■  J  1 

(23) 

V'nl  =  2  do,E^ 

(24) 

vfl  =  4 

(25) 

Vfl  =  2d&e0{Ep  +  2E^EW) 

(26) 

rtfl  -  idhoE^EW* 

(27) 

vil  =  4  d4t0(E^EW'  +  E'3)£f  *) 

(28) 

-  4dfie0(45)^0)* + q3)q2)* + Bi1)£j4)*) 

(29) 

Using  equations  19  and  24,  and  substituting  in  18,  the  first  term  on  the  right 
hand  side  of  eqn.  20  can  be  determined.  Using  this  expression  for  dE) 1  and 
eqn.  19  in  eqn.  27,  the  dr term  in  eqn.  21  can  be  evaluated  from  eqn. 
17.  With  this  d?Ep}  and  the  Es  from  eqn.  19  substituted  in  eqn.  25,  equation 
17  can  be  used  for  the  evaluation  of  the  dsE*^  term  in  Eqn.  20,  and  so  on. 
Before  embarking  on  the  explicit  determination  of  each  of  these  terms,  in  the 
next  section  several  parameters  are  defined,  and  a  new  notation  is  introduced 
which  allows  ail  the  spatial  integrations  to  be  performed  in  dimensionless  units. 


3.1  Notation 

The  parameter  £p  =  i/bp  where  l  denotes  the  length  of  the  nonlinear  medium 
is  used  to  characterize  the  amount  of  focusing  of  the  pump  beam.  The  phase 
mismatch  between  the  nonlinear  polarization  beam  and  the  generated  SHG  wave 
is  denoted  by 

Ak  =  2kp-ks.  (30) 

Similarly,  the  mismatch  in  linear  absorption  between  the  nonlinear  polarization 
and  the  SHG  beam  is  proportional  to 

AasOp-0,/2.  (31) 

In  addition,  the  following  parameters  are  defined: 

p  =  x/wop  q  =  y/wop  t  =  zfl 
Pi  ~  x  j wop  q\  =  yf l h  =  z*  jl 


5 


(32) 

(33) 


Pa  =  Xa/Wop  Wl  -  pl/w 0p 


(34) 


where  W\  is  a  dimensionless  parameter  which  provides  a  measure  o£  the  amount 
of  walk-off  of  the  SHG  wave  as  if  propagates  through  the  crystal.  Also,  defining 
the  variables  and  dp(t)  as 


Qp  —  1 

di(t)  =  ap  +  2i£pt 


(35) 


the  incident  pump  field  can  be  expressed  as: 

£<0)(r)  _  EP°  C-(V2  +Q!)/d1{t) 

p 


(36) 

(37) 


4  The  expressions  for  the  electric  field  terms 

4.1  The  first  order  SHG  field  term 

Using  Eqns.  18,19  and  24,  the  first  order  signal  field  term  dE^ 1 1  is  obtained 
(after  integrating  over  the  transverse  coordinates) : 


d£W(r)  =  *£w>/il)(r) 


where 


and 


2mdlETi$ 

K  =  —  - 


Jo 


A  £  n  g 


e(iAfc- Act)  Hi  - (P* 2 +</a ) 


Here  pa  —p  -  W\(t  —  )*  The  function  ^(Mi)  is  defined  as 

d^it,  ti)  =  di(ti)  +  2 ik£p(t  -  ti), 


(38) 

(39) 

(40) 


(41) 


where  k  =  2  kp/ha. 

4,2  The  second  order  pump  field  term 

Using  Eqns  38,19,  27  and  17,  the  second  order  pump  field  term  <pEP2}  can  be 
obtained  (after  integrating  over  the  transverse  coordinates): 

^(r)  =  -^\K\2f^(r)  (42) 
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where 


^  c(ifeP-ap/2)(*  J  ^^e(-iAk-ot9/2)ltl 

X  jU  dt 

g  -  (p2  +  g2  )02  +2p62  +c^ 

(43) 

dl(i4)d8(t,tl,  t2>^3j  ^4) 

with 

a2(t,ii,t2)-  d4(Mljt2y 

(44) 

(45) 

and 

(46) 

Other  functions  used  above  are  defined  as 

d2a(f,  U)  =  di(ti)*  4-  2i^(£  -  ti) 

(47) 

dz(tut2)  =  is}  +  2di(£j)* 

(48) 

d4(t,ti,  £2)  =  di(£i)*d2(ti,t2)  +  2i£p{£  -  fi)d3(ii,  ta) 

(49) 

4.3  The  third  order  SHG  field  term 

Using  Eqns.  19,  42,  25  and  18,  the  third  order  SHG  field  term  d3E[3']  can  be 
obtained  (after  integrating  over  the  transverse  coordinates): 


where 


where 


d3 £^3) (r )  =  -(— }2rc  ]  k  [2  /a(3)(r) 


^(3)(r)  _  e(ifc.-a./2)K  f*  dtie(i&k-ba)U, 

8  Jo 

x  F  dt2e(--^k-a-/^  r  dt3eliAk-Aa)lt> 

J  o  J  o 

^  -  (p*2  +S2  )03  +  2rja  +C3 


di  (^3  )ds(tj  fii  £2) 


ji,  t2)  ^3) 


(50) 


(51) 

(52) 


7 


(53) 


and 


with 


and 


=  2lFi  (t'2  —  ta) 
C3(*,ti»*a)  =  ^i2(^  _  ^s)2 


df>{tytut2,hy 

gi(Mufrz) 


t\y  £2)^6  (tj  tu$2->  £3) 


(54) 


<^5  (^1 5  ^2  ?  ta)  =  ^4  1 3  ^2 1  £3)  +  di(£i)d3(£2,  £3)  (55) 

dfl{Mi»*2its)  =  di(ti)d4{ti,t2,t3)  +ik%p(t  -  ti)d5(^i,t2)  (56) 


—  4ik£p(t  -  ti)di(ti)di(t2)*2  -  2d2a(ti,  t2)d$(ti  ti ,  ta,  £3)  (57) 


4 A  The  fourth  order  pump  field  terms 

Using  Eqns.  38,  42 ,  50  and  19  in  28  and  17,  the  fourth  order  pump  field  term 
d4E^]  is  obtained  as  a  sum  of  two  terms  and  d^E^V'. 

*E$(t)  =  2(“)2  |  k  |4  E^flfir)  (58) 

UP 

rf444)(r)  =  (Ji)2|«|4£po42){r)  (59) 

np 

The  expressions  for  and  are  given  in  the  next  two  subsections. 


4.4.1 


C(r)  ^  e{ik 


where 


r-<Xy/m  f  rftiC(-iAfc-a./2},ti  f  1  ^(iAk-A^ 

J  o  Jo 

f2  dhe(-iAk~a''2)ltl  f3  dt^iAk~Aa)lti 

Jo  Jo 

g  —  (pa  +  93)^41  +2p^Jl  +  £<tl 

di  (^4  )dg  (f,  ij. ,  t-2 » tsi  ^4 ) 

,  ,  ,  ,  .  _ 

a4i(t,ti,*a,t8. *4)-  dg(Ml)(3,t3ii4)’ 

_  ...  di(*i)*ga(ti,*2,t3,t4) 


641  iji  ti,  ia,  £3,  (4)  —  W'*i 


^3  (£1  f  1  j  ^2 .  ^3 1  ^4 ) 


and 


(60) 

(61) 

(62) 


C4l(Ml,i2,f3,i4)  =  W12{ff3(tl,i2,f3,*4)  +  ffl{*itl.*2.f3)t4)}  (63) 
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with 


92{tuhi  £3*  *4)  =  (£1  —  £2)^5  (£21  £3*  h)  -  %{h  -  £4)^1  (£3)*^  (£2)  (^4) 


^3(^1 1  £2*  hi  £4) 
94(^7  ^li  hi  hi  *4) 


ff3n.(£li  £2?  ^3  1  ^4) 
<?3d(ii  j  hi h,  £4) 

94n{t,  hi  hi  hi  £4) 

94d{ti  ti,h,  £3,  £4) 


(65) 

(66) 


9zn(tu  hi  hi  £4)  —  {h  —  *4)2tfi(ii»  hi  £3*  £4) 


-(£1  -  t2)2^4{hi  hi  i4)ds(£2>£3i  £4)  (67) 

93d{tu  hi  hi  h)  —  d^lhi  hi  h)d^{ti^  hi  hi  h)  (68) 

$4n(t» tiitzi taiM  =  2 i£P(t  -  t^d^tiY g2{tut2it2,h)2  (69) 

g4d(thhi  hi  h)  —  ^6 (iii  hi  hi  £4)^8  (£1  £i*  £2*  £3*  £4)-  (70) 

The  functions  di  and  de  are  given  by 

di (thhihih)  =  d$ (i  1 1  hi  hi  h )  +  di (t  1 ) * d^ih j  £3  >  £4 )  (71 ) 

and 


4 (£1 ,  £2 1  £3 1  £4 )  =  4  (£1  )*4 (£ii  £2 » £3 1  £4 )  +  2 i^P{t  —  £1  )df  (£1 ,  £2 1  £3 » £4)  (72) 


4.4.2 


where 


and 


a^it  fdt ,e< 

Jo 

f  l  dhe{iAk~ajm*  f*  dtAe(^ 

Jo  Jo 


=  e(lfcp 
rt  1 


-iAfc-ac/2}tt| 

^3 


ft\ 

/  dt2JiAk 

Jo 


—  AfiEjftj 


— tAfc— Aa)ft4 


p  — (p  +q  )<l42+2p£)42  +  C42 

dj:(£4)dg(£,  £i?  £21  hi  £4) 


a42(£,£l,t2,£3,£4) 


da(£i, £2?  £31  £4) 
dio(£i  £i?  £sj  £3?  £4) 


&42(£j  £1?  £21  hi  £4)  — 


di{ti)*g^{ti,hihih) 
dio(£t  £ii  £2*  hi  £4) 


C42(£,£i,  £2?  £3,  £4)  =  W/,i2^e(£ii  £21  £31  £4) 


(73) 


(74) 

(75) 

(76) 
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with 


h,U)  =  2(ti  -  £2)d4(£i,£a,£4)*  +  2(t3  -  U)d2i(tiit3)*d2(tut2)  (77) 


0e(£»  £i>  £27  ^3t  £4)  = 


2igp(t  -  tl)95(tl,t2,t3rt4)2 
f2)d|(ti,  £3,  £4)*rfio(£,  fli  ^2t  ^3 j  £4) 
2(ti  “  £2)^  2(£3  -  £4)^21  (ill  £3)* 


(78) 


^2(^17  £2)  <£j(£ii  £35  £4)* 

The  functions  dg  and  d10  are  given  by 

^(£17  £21  £37  £4)  =  2^4  (ti ,  £35  £4  )#  T  da  (£  1  >  £2)^3  (£37  £4  }*  (79) 

and 

f£io(£i  £  1. 7  £2,  £37  £4)  =  d2(£i,  £2)^4 (£ii  £3*  £4)*  +  2i£p(£  -  £i)dg(£|,  £27  £37  £4)  (80) 


4.5  The  fifth  order  SHG  field  term 


Using  Eqns  42,  58 ,  59  and  19  in  26  and  18,  the  fifth  order  pump  field  term 
is  obtained  as  a  sum  of  three  terms  d5£^j\  d5E^  and  #1?^: 


=  (r1)2^  I  *  I4  EpoflV  W 


(5)/ 


d&EH}{r)  =  4(^-)2k.  |  k  [4  Epoflzir) 


and 


d5E%( r)  =  2(^)2«  |  k  |4  Ep afi»{r) 


(5)/ 


(81) 

(82) 

(83) 


where  the  expressions  for  the  terms  /s(\  /^ '  and  are  given  below: 


4.5.1  /<? 


fal\T)  = 


=  e{ik9~aa/2)tt  f  ^ e(i&k-Aa}ltt 
JO 

f  '  rft3e(iAfc“A")ii3 

Jo 

[  1  dUe(~iAk-a’/2)lu  f  1  dt5e{i 
Jo  Jo 


x  I  1  dt2e{~iAk~a‘l‘1)iti  r  2 
Jo 


(tAfc  — 


e  ~{pa  2+q2  )a$i  +2pa  +c3i 

(£3)^1  (£5)^12  (£j  £1 7  £2!  £37  £47  £5) 


(84) 
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where 


(85) 


ae,  (t,tut2,t3>t4)=  MM***) 


and 


di2{t ,  ^1,^2, ^3, *4*  is) 

t.,(M„fa.t„<,).iy,/’(‘1;‘,,'¥irM, 

“12^5  Mr  Mi  Mr  Mr  Ml 


_  tji r  2  rtfs(£]  j  *4?  4)2^12a(*i  ill  ^2i  *4r  is) 


£si  (i,  i 1 1  ^2 ,  ts r  U )  —  W\  { 


^12  ($1  ^1  r  ^2  ]  ^3 r  ^4  ]  is) 

£2,  ^4 r  is)} 


with 


97(h,  t2>  ia,  ^4,  is)  —  2(ia  -  i3)^i(i2)*d4(ii,  i4i  is) 

4-  2(^4  —  ts ) c?i  (^4 ) * ^4(^1 5  f2i  ia) 

9&(tu  *2t  is*  i4i  is)  =  %(h  _  is )  ,  , !  "  , — h  ^(£4  ~  £5)  -  - •  77  •  / ''  ,  -. 

2  ^21  (ti ,  t2)  f  .2  ^21  (i  l?  i4.) 


gg(tl,t2,  i3i  i4r  is)  =  2(^2  -  is) 


+  2(ti  “i5r 


^4(il 1 i4r i&) 


^4(iir  i2i  is 
The  functions  dn  and  dn  are  given  by 

dn  (iii  £21 4,  i4)  =  dajia,  ia)d4{ii,  tir  h)  +  ^3(^4  j  is)d4(ii,  is) 


(86) 


(87) 

(88) 

(89) 

(90) 

(91) 


and 


di2(i»  iii  ^2r  i3j  ^4)  =  ^(ii^isiisJ^Ciii^jis) 

+  ik£p(t  -  fi)dii(iiT£2lt3ri4,is)  (92) 


where 


(5) 


4,5.2 


/^(r)  =  e(ifc*_£“'/2),f  [  dt 

J  0 

x  f  *  dt2e{~i&k~°-mtl  [  ^  dt3e(i 

J  0  Jo 

x  [  '  dt4e<--iAk-a~mu*  (U  dt5e(i*k 

Ja  J  0 


1  Ah— Aat)lt$ 


”4a)lt  5 
g — (Pa  3  +9  3 ) «a  2 + 2pw  &2  +  Cfia 

di(£a)di(ta)di2{i,  iir i2r  is?  i<ii  is) 
dia(ii, (21 iai *4* is) 


^5a(irili  i2r  ^3 1  £4)  — 


di4{i,  ii,  i2i  i3,  i4i  is)  * 


(93) 

(94) 
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(95) 


and 


^52(^7  ti,  fa,  fa,  fa)  =  W\ 


di(ti)di(fa)*  g2{fa,fa,t4,fa) 

dw(t,  1 1 ,  ta,  t3,  £4,  £5} 


cs2(£,  fa,  fa,  £3,  £4)  =  ^i2{^3(£2i  £3,  £4?  fa)  +  g^i,  fa,  fa,  £41  £5) 

+  i*€P(t-ti)di(ti)di(t3raAi(ta.*3,t4.Ma] 

^s(£ii  £21  fa,  £4*  £s)^m(£i  £i?  £21  £31  £4?  £5)- 
The  functions  and  di4  are  given  by 

disifayfa,  fa 1  £4,  £5)  ^  (£1)^7  (£2  ,£3,  £41  £5)  +  tf g  C  £  1 1  £  2 :  £31  £4^  £5) 

and 


(96) 

(97) 


du{£i  £n  £21  £31  £4?  £5)  =  (£1)^3  (£ii  £2,  £3 1  fa,  £5) 

4-  ik£p(t  —  £i)dj3{£i,£21  £31  £4*  £5)  (98) 


4*5.3  /., 


(5) 


/i?W-«Ci*'““'/a)l<  /  dfae{i^k 

Jo 

x  f  '  d£2e(“iAfe"Qa/2)ifa  [  2  dfaemk 
Jo  Jo 

X  f2  dfaeiiAk~«‘/2)lt*  fU  dfae(- 
J  0  Jo 


-&a)Ui 

-  A»)itg 


where 


iAk  —  £ct)ltz 


e~(p*  )aaa+2pal>a3+co3 

d\ (£3)^1  (t5)<fi2(i,  fa,  fa,  fa,  fa,  fa) 

{,  ,  ,  ,  ,  di5(fa,  £2l  £3,  £4*  £5) 

a53l£j  Hi  fa j  fa,  fa)  — 


and 


tfie(£>  £1,  £2  t  £3 1  £4 ,  £5 ) 

l  tx  ,  j.  ,  ±  \  -  117  di{fa)9${fa,faifa,fa) 

^53^,  fa,  fa,  fa,  fa)  =  Wi  2  {<?6  (£l »  £21  £3 1  £4?  £5) 

^  ^  ( £  £i)tii(£L)  £t3  { £  2  *  £35  £4^  £  s ) 


(99) 

(100) 

(101) 

-}  (102) 


dio(ti,fa,fa,  fa,fa)dw(t,  £i,  £2,  £3?  £4,  £5)^ 

The  functions  di5  and  d|6  are  given  by 

f£is(£ii  £2^  fa,  fa,  £5)  =  {£  1  )<^9 (£2 j  £3*  £43  £5)  +  t£i_o(£;t,  £2,  £3,  £4,  fa)  (103) 

and 

dietMii  £2?  £3,  £47  £5)  —  ^i{£i)'^io(£ii  £27  £33  £47  £5) 

+  ik£p(t  —  £i)di5(£j,  £2,  £3 7  £4?  £5)  (104) 
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5  Conclusions 


Expressions  for  the  pump  and  second  harmonic  fields  are  presented  as  series  of 
terms  in  increasing  powers  of  the  nonlinear  optical  coefficient.  The  effects  of 
beam  diffraction,  linear  absorption  and  beam  walk-off  arc  included.  An  advan¬ 
tage  of  the  method  presented  here  over  the  widely  used  numerical  method  may 
lie  in  the  increase  in  the  computation  speed,  especially  in  special  cases.  For 
example,  in  the  case  of  a  collimated  pump  beam,  i.e,,  for  1,  the  single, 

double  and  triple  integrals  in  the  expressions  for  the  fields  reduce  to  analytical 
expressions  and  the  quadruple  and  quintuple  integrals  reduce  to  much  simpler 
integrals  which  can  be  rapidly  evaluated  using  standard  computers.  Also,  the 
method  can  be  extended  to  the  case  of  pump  beams  of  arbitrary  shape. 
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